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Coherent production of lepton pair in muon scattering by nucleus is considered. 
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O ; With the help of different approaches (equivalent photon approximation, direct nn- 

merical and analytical calculations) we obtain total cross section of this reaction 
and distributions over invariant masses of final leptons. Special attention is paid 
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to influence of final lepton identity on total and differential cross sections. We also 
study the role of coherence condition. Total cross sections of charmonia production 



in muon interaction with electromagnetic field of nucleus are also presented. 
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(N ! I. INTRODUCTION 

In our article we consider total and differential cross sections of lepton pair production 
in muon scattering by nucleus 

O 



fi + M -> (1) 

This subject is discovered thoroughly already (see, for example, , but there are 

some questions, that deserve more detailed investigation. 

One of this questions is the influence of final lepton identity on total and differential 
cross sections. In mentioned above works it was shown, that total cross sections are changed 
slightly when this effect is taken into account. Since final lepton identity strongly complicates 
calculations (for example, the number of Feynman diagrams is doubled), in the following it 
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was often neglected. Actually, the following process was considered 

with mass of £-lepton equal to that of muon. In our article we will show, that, though final 
lepton identity does not change total cross section significantly, invariant mass distributions 
are changed dramatically. This difference is important for calculation of electromagnetic 
background to charmonia production in muon scattering by nucleus. 

The second question is the influence of coherent condition (that is the requirement that 
the nucleus stays undestroyed after muon scattering). This condition can be described by 
nucleus electromagnetic forfactor, but there is not enough information about its explicit 
form. In the framework of equivalent photon approximation (EPA) method, that will be 
used in our paper, there is a simple way to check this condition. One only needs to take a 
suitable distribution function of equivalent photons. 

In second and third sections of our article we calculate total and differential cross sections 
of the process §B) using three different approaches - EPA, analytic expression, given in fl, 
and direct calculation. In fourth section the influence of coherence condition is studied. 
Finally we consider charmonia production in muon interaction with electromagnetic field of 
the nucleus. 



II. TOTAL CROSS SECTION 



Leading order diagrams for process 

H + {k)M(p) - ^ + (k^(p r )£ + (k i+ )£-(k e -) (2) 

are shown in figures (TJ [2j If there are identical final particles (i.e. £ = fi) one should also 
add diagrams with permutation of £ + and /i + lines. 

First of all, it is clear that diagram fig. [lb does not give significant contribution to the 
cross section, since corresponding amplitude is suppressed by nucleus mass Mj^. It is also 
clear, that amplitudes of other diagrams are large at small virtuality of the photon emitted 
by nucleus. So when calculating the cross section of this reaction one can use Equivalent 
Photon Approximation (EPA) . In the framework of this method the cross section 
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FIG. 2: ^77 vertex 

of the process (j2j) is expressed through the cross section of lepton pair production in muon- 
photon interaction: 

V+ikMq) - /i+(VK + (^+K~(V)- (3) 

The diagrams describing this process can be obtained from diagrams fig. []Ji,b after removing 
nucleus current. 

According to equivalent photon approximation cross sections of reactions ([2]) and ([3]) are 
connected by the relation 



J a^(m)n(m)drh, (4) 



where m is the invariant mass of final leptons (m 2 = (k e - + k^+ + k^+) 2 ) and n{m) is the 
distribution function of equivalent photons in nucleus. In logarithmic approximation this 
function is 1 -* 

, „ N , _ 4Z 2 a ( m s \ dm 
n(m)dm = In I — — — I — — , (5) 



7r V Mj^f m 2 J m 



1 ) The difference between this expression and function given in [7|] is explained by the variable transforma- 
tion. We use the invariant mass m instead of virtual photon energy 
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FIG. 3: Total cross section of fij — * 3/i reaction as a function of rh 

where a is fine structure constant, Z is nucleus charge number, s = (k + p) 2 , m and M^y are 
muon and nucleus masses respectivelly. The integration domain in (jlj) is determined from 
/ry — > 3/i threshold and the condition of positivity of distribution (j5j): 

(3m) < rh < —s. 
K ' M 

In fig. [3] we show yU7 — > 3/i cross section calculated with numerical integration 2 ) as a 
function of rh. Using this cross section and formulae (jlj), (jSJ) we have computed total cross 
section of muon pair production in muon scattering by lithium (that is the nucleus with 
charge and mass numbers equal to Z = 3 and A = 7). It is interesting to mention, that for 
fixed muon energy in laboratory frame (that is the rest frame of initial nucleus) this cross 
section does not depend on nucleus mass. So it is convenient to introduce a dimensionless 
variable 

s ^ 

2 - ) In our calculations we have used two methods: COMPHEP 0] and method described in [9]. The difference 
between them is that in the first case squared matrix element is computed analytically, while in the second 
method matrix element computation is done numerically and the resulting value is squared. It is clear, 
that for large number of diagrams the latter approach requires much less numerical calculations and, 
hence, leads to smaller errors. The results of these methods coincide. 
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where 7 is the lorentz factor of initial muon in laboratory frame. 

This cross section should be compared with the results of analytic calculation of /jj\f — > 
fiAffifi cross section. In work 5j the reaction ep — > epee was considered. In the logarithmic 
approximation the cross section of this reaction is 

a{ep epee) = fr 8 - AL 2 + BL + C), (6) 

where 



L = In -— - , (7) 

\m e M p ) 1 ; 

and coefficients A, B, C are equal to 

A « 6.36, B « 2.6, C « 40. 

For large s these coefficients do not depend strongly on particle masses, so for our case we 
should only change lepton and nucleus masses in the above expressions. We would like to 
mention, that, analogously to EPA, for fixed muon energy this expression does not depend 
on nucleus mass. 

Another way of calculation of ijM — > fiAffifi cross section is direct numerical calculation. 
In this case we have observed the dependence of cross section on nucleus mass. If Mv and 
s are increased (while muon 7-factor is kept constant) total cross section of (j2J) decreases. 
This decrease is caused by power corrections (~ Mjj/s) that are neglected in EPA and ([6]). 

In fig. H]we show fxAf — > 3/i/V total cross section as a function of x. It can be seen, that 
the results of direct calculation (dots) and formula ([6]) (dashed line) are close to each other 
(the difference for low x is explained by power corrections), while EPA cross section (solid 
line) is smaller. This discrepancy however, can be removed if we multiply the argument 
of the logarithm in (jSJ) by ~ 1.4 (dotted line in fig. |4|). Since formula (jSJ) was obtained 
in logarithmic approximation, this argument is defined up to order one factor, so such 
transformation is allowed. 

Let us discuss now the influence of final lepton identity on total cross sections of reactions 
(J2J) and (jSJ). The amplitude of the process /i + 7 — > is obtained from diagrams shown 

in figOJi,b (after removing nucleus current and anti-symmetrizing identical final lepton legs) 
and can be written as 
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FIG. 4: Cross section of //Li — ► 3/zLi reaction. Solid and dotted lines are EPA results with 
distribution functions ([5]) and corrected distribution function (see text for details); dashed line is 
the result of formula ©, dots — direct calculation. 

where A^[ i s the amplitude of process ([3]) without final lepton identity, and symbols [i and 
I designate momenta and helicity variables of and £ + -leptons respectively. Total cross 
section is equal to 

a? = yd^\An 2 = yd^{\A^J)\ 2 +\A^(e,fi)\ 2 }- 

- lfd*{^MA&%i,ri + ^*M^(e,iM)}, (8) 

where integration is performed over the loretz-invariant phase space of final particles (aver- 
aging over their polarizations is also implied), and factor 1/2 is caused by identical leptons in 
final state. Since integration domain is symmetric under permutation /i + <-» £ + , expression 
(jSJ) can be rewritten as 

of = J d*\A&(ji,£)\ 2 - Re/d$ {A&Qi,£)Agr(l,ri} = 

= a£l + 5a^. (9) 

First term here is the cross section with final lepton identity neglected, and second term 
describes the corrections caused by this effect. It is easy to understand, that for large s this 
correction is small. In the region of low invariant mass m^+^ = J {k(+ + kg-) 2 the virtuality 
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FIG. 5: The ratio of fxj — > 3/x cross sections, calculated with and without final lepton identity 

of show in fig. [Tb photon is small and the amplitude A^{(£, /i) is enhanced. So in this region 
both factors in the first term of expression are large. This is not the fact for interference 
contribution (second term in 0) since there is no enhancement for A^{(1, fi) in this region. 

We would like to stress, that presented above arguments do not depend on the dynamics 
of specific process. The only requirements is symmetry of phase space integration domain 
under permutation of final particles and an enhancement of the amplitude on the border of 
this domain. Near the threshold the interference term will be comparable with uf 7 and final 
lepton identity will play significant role. This is because for low m total integration domain 
is comparable with the width of the region, where amplitude is enhanced. 

In figure [5] we show the ratio of fij — > 3/i cross section with final lepton identity taken 
into account and neglected. This figure agrees well with out expectations. Since in EPA 
cross sections of reactions /ry — > 3/x and fj,J\f — > 3/zA/" are related by formula (@J, for the 
latter process the behavior of such ratio will be the same. 



III. INVARIANT MASSDISTRIBUTIONS 

Let us now discuss distributions over invariant masses 

m 2 = (k e - + + k i+ ) 2 , 
m l+£+ = (V- + h+) 2 , 



8 




(k^+ + k e -) 2 , 
(k e+ + k e -) 2 . 



In the case of identical leptons last two variables coincide. 

In the following we will consider the scattering of 160 GeV muon by lithium. The dis- 
tribution over mass rh follows directly from eq. (j4j). In figure [6] we show this distribution 
calculated with EPA (soled line) and direct calculation (dots). From previous section (see 
figure E]) it is clear, that final lepton identity has little effect on this distribution. 

Distribution of muon-nucleus cross section over invariant mass of /x + £ + -pair in EPA is 
calculated with formula similar to (jl]): 



This distribution and the ratio of differential cross sections with and without final lepton 
identity are show in figure [71 It is clear, that in this case the influence of identity is also 
small. The reason is clear — the integration domain is symmetric under the permutation of 
/i + and £ + momenta. 
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FIG. 6: Distribution of /iLi — > /iLi/ifi cross section over invariant mass of three final leptons 
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FIG. 7: (a) — distribution of juLi — * /iLi/i/i cross section over m^+£+; (b) — the ratio of distribu- 
tions with final lepton identity neglected and taken into account 

If we consider the distribution over invariant mass of leptons with opposite charges (i.e. 
the situation is completely different. When final lepton identity is neglected, we 
actually have two different variables — mg+i- and m^+i-. If identity is taken into account, 
theses variables coincide. In figure [8] we show distributions over these three masses. It can 
be clearly seen (see fig. (U where the ratio of distribution with and without final lepton 
identity is shown), that final lepton identity changes these distributions dramatically. If we 
measure, for example, the distribution over pair mass in the region of Mj/^ with an 

instrumental error A = 50 MeV, we get 

da^ 



A- 



dm^i 



1.5 pb (10) 



when final lepton identity is neglected and 

. da^ 



dm^i 



92pb. (11) 



- m j/ii> 



if it is taken into account. This difference is crucial in determining the electromagnetic 
background to J/ip production in muon scattering by nucleus. 



IV. FORM-FACTORS AND COHERENCE CONDITION 

Up to now we treated nucleus as a point-like charged particle. It is clear, that this 
approximation has little to do with real world, since nucleus can easily be destroyed in 
muon scattering. 



10 



da nb 

dm 4. _ GeV 




ij,- j GeV 



FIG. 8: Distributions of //Li — > /iLi/U/x cross section over invariant mass of leptons with opposite 
charge. Solid line — distribution with final lepton identity taken into account; dashed and dotted 
lines — distributions over m^+i- and mg+£- ; dots — results of direct calculation 
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FIG. 9: Ratio of /xLi — > /iLi/i/i cross section distributions with final lepton identity taken into 
account and neglected 

One can go beyond this approximation by introducing electromagnetic form factor of the 
nucleus. It is well known (see, for example, [h]] ), that for proton this form factor is 

F{q2) = ' Qo ~ 0.7 GeV 2 . 

For nucleus formfactor is unknown, but we can describe it qualitatively by imposing a 
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restriction on a squared transferred momentum — q 2 < Qq. In our calculations we have used 
values Ql = 0.1 GeV 2 and Q 2 = 0.05 GeV 2 . 

In the framework of equivalent photon approximation it is more convenient to use another 
approach p . The dependence on impact parameter b is introduced in distribution function 
n: 

where u is photon energy and i£o,i are modified Bessel functions. For nucleus to remain 
undestroyed the following condition should should be valid: b > R mm . Here R min rs 1.1A 1 ' 3 
fm is a typical radius of the nucleus. After integration over impact parameter with this 
restriction we get 



, s 2Z 2 a 
n K (u) ■■ 



TTUJ 



£jr o (0tfi(O-y (k 2 (O-k 2 (O 



(12) 



where £ = uR m i n / ^(3. 

In figure [TO] we show distribution of fiLi — > fiLififi cross section (with final lepton 

identity taken into account), obtained with EPA (distribution function f)12p was used) and 
the results of direct calculations for different values of a cutoff parameter. 



V. ijM [i Xc M 

Let us now proceed to charmonium production in muon scattering by nucleus. In par- 
ticular, we would like to consider x c -meson production via the interaction of muon with 
equivalent photons (corresponding diagrams are shown in figure HTb) . In comparison with 
gluon reaction (fig. fTTb) the cross section of this process is suppressed by electromagnetic 
constant. In the gluon case, on the other hand, color conservation requires emission of an 
additional gluon. It is therefore straightforward to get the following estimate for the ratio 
of electromagnetic and strong cross sections: 

<r(q9 -> QXcJ) a 2 

In paper [ll| one can find the expressions for differential cross sections of charmonium 
production in quark-gluon interaction. It is evident, that in our case we should remove from 
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FIG. 10: Distribution of fihi — > fihi^ifi cross section over invariant mass of (/i + /U~)-pair. Solid 
line is EPA result with distribution function eqrefeqmK, circles and squares are represent direct 
calculation with — q 2 < 0.1 GeV 2 and — q 2 < 0.05 GeV 2 respectively, dashed line is EPA in point-like 
approximation . 




FIG. 11: Diagrams of Xc- m eson production in (a) electromagnetic and (b) strong processes. 

these expressions all color factors (including factors caused by color degrees of freedom of 

initial particles) and change a s to a. As a result we get 

dajpg -> fix) nAT « 3 dcrjqg -> qx) n ^ 
dq~ 2 " 27Vc ^| dq~ 2 • [l6) 

There is another difference between these two reactions. In our case, contrary to scattering 
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FIG. 12: Solid line — J = 2, dashed line — J = 0, dotted line — J = 1. 

on light quark, we cannot neglect muon mass. The equivalent photon approximation, used 
used in our paper, give results with logarithmic accuracy, so there is no need to leave muon 
mass in the differential cross section (1T3"|) . The inffuence of this mass on the integration 
domain is, however, crucial. The lower bound on the squared tranfered momentum is 



where M x is x<r m eson mass (we neglect the difference between Xco, Xci and Xc2 masses). In 
the limit = we have t m i n = 0. Cross section (TlBl) diverges logarithmicaly in this region, 
so muon mass should be taken into account. 

In figure [T2] the dependence of the /i7 — > /i% cross section on /i7 energy is shown. After 
integrating with distribution function of equivalent photons we obtain the cross section of 
the reaction /zLi — > fiLixj: 



o.oi 



Miml 



(14) 



t 



mm 




cr(/iLi — > /iLix c0 ) = 14 fb, 



(r(yuLi — ► /iLix c i) = 0.5 fb, 



cr(/iLi — > /iLix c2 ) = 19 fb. 
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Produced x-mesons decay subsequently into experimentally observed lepton pair (xc — > 
— > yU + /i _ 7). If we neglect permutation background, than the cross section of this 
process is 

<r(/iLi — > fiLix — > fiLijfi + fi~) = a{ji\A — > /iLix)Br(x — > ■yJ/ip)Br(J/ip — > n + 

For different values of Xc-meson spin we have 

cr(/xLi — > /LtLiXco - ► /iLi7/i + /i~) = 8 • 10~ 3 fb, 
cr(/iLi — ► fiLiXd /iLi7yU + /i") = 11 ■ 10~ 3 fb, 
cr(/iLi — ► fiLiXc2 — * /iLi7yU + /i~) = 0.22 fb. 

It is clear, that these cross sections are smaller, than the cross section of the non-resonant 
process (JTJ). 



VI. CONCLUSION 



In our article we consider lepton pair production in muon scattering by nucleus. This 
process was already studied in details in earlier works, but some questions are still opened. 

One of such questions is the influence of final lepton identity on total cross section of the 
process 



and distributions over invariant mass of yU + /i~-pair. In earlier works it was shown that 
this identity changes total cross section slightly. We confirm this result. Invariant mass 
distributions, on the other hand, change significantly when final lepton identity is taken into 
account. For example, if muon with energy 160 GeV is scattered by litium, than in the 
region m M+M - 



Mj/j/j we have 



da(/Aji — ► /iLiyU + /i ) 



c?m At+jU - 



U 



m=M 



J/ip 



nb 
'GeV 



for identical leptons and 



da(fiLi — > fiLifi + fi 



dm 



0.03 



m=M 



J/4, 



nb 
GeV 
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when this effect is neglected. This difference is significant when estimating background to 
J/^-meson production in lepton- nucleus scattering. 

In the final section we used Weitzaker- Williams method to calculated the cross section of 
Xc-meson production in electromagnetic interaction of muon with nucleus and subsequent 
decay x ~^ ~~ 7/i + /i~- Resulting cross sections are significantly smaller, than the 

cross sections of the nonresonant processes. So, one can state that charmonium mesons in 
lepton-nucleus interaction are produced by strong interaction. 
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